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Abstract 

Considering the standard abelian sandpile model in one dimen- 
sion, we construct an infinite volume Markov process corresponding 
to its thermodynamic (infinite volume) limit. The main difficulty we 
overcome is the strong non-locality of the dynamics. However, using 
similar ideas as in recent extensions of the standard Gibbs formalism 
for lattice spin systems, we can identify a set of 'good' configurations 
on which the dynamics is effectively local. We prove that every con- 
figuration converges in a finite time to the unique invariant measure. 
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1 Introduction. 



The standard abelian sandpile model was introduced in 1988, |l|]. This model 
lives on a finite graph (e.g. a square on Z 2 ) in which to each vertex is 
associated a height-variable (representing the height of a sand column at 
that site or, the averaged difference in sandlevels with neighboring sites). A 
vertex is picked at random and its height is increased by one. If the height 
variable exceeds some critical value, then the vertex becomes "unstable" and 
"topples", i.e. gives an equal portion of its sand to each of its neighbours 
(adjacent vertices) which in their turn, can become "unstable" and "topple", 
and so on, until every vertex has again a subcritical height-value. Some 
vertices are connected to a sink (the boundary) where sand disappears. An 
unstable vertex thus creates an "avalanche" which can cause the toppling of 
many vertices around it. The number of vertices involved in or affected by 
one avalanche can be arbitrarily large (dependent on the configuration). 

Over the last decade, sandpile processes and various versions thereof have 
been intensively studied as simple lattice models exhibiting the phenomenon 
of self-organized criticality. While it is still a subject of debate what is 
the precise nature of the phenomenon and under exactly what conditions it 
obtains, sandpiles have become widely appreciated as simple threshold dy- 
namics showing scale invariance of various response functions with (critical) 
exponents of an astounding universality. Not surprisingly, a vast amount of 
computer simulations have been performed involving sandpile processes and 
their acclaimed self-organized criticality. Yet, both controlled laboratory or 
real nature experiments and a rigorous mathematical analysis have remained 
rather limited. We refer to |L3| and the references therein for a number of ex- 
periments validating some of the universal aspects of the sandpile paradigm. 
In this paper we do not address questions of self-organized criticality nor do 
we refer to the most interesting physical realizations of sandpile-like models. 
Rather we turn to the simplest mathematical questions one can ask from the 
point of view of the theory of interacting particle systems for the simplest 
one-dimensional abelian sandpile process. 

For the abelian sandpiles the first mathematical results have been ob- 
tained by D. Dhar, see || [| f|] and the review paper 0. As further discussed 
also in |Tj| , the main question there was the characterization of the invariant 
measure and the recurrent configurations of the finite Markov chain defined 
in the sandpile process. In one dimension, the answers are rather simple. 
A detailed study can be found in [TH|. Our questions concern the so called 



thermodynamic limit, i.e. how to construct a bona fide infinite volume limit 
dynamics and how to characterize its main properties. In [O] , Priezzhev was 



able to compute exactly the single-site marginal of the invariant measure of 
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the sandpile process on a square in the limit as the square tends to cover the 
whole plane. Various other mathematical results have been obtained but, for 
all we know, the problem of the thermodynamic limit has remained widely 
open so far. In the present paper we give a complete solution to this problem 
for the standard abelian sandpile process in one dimension. 

It should be realized that this one-dimensional model has as such rather 
limited interest from the physics point of view. Yet, some of the mathematical 
questions pertaining to the thermodynamic limit remain very non-trivial and, 
arguably have answers which are expected to be not so different in higher 
dimensions. It is therefore unavoidable to start in one dimension and clarify 
the situation there. The main mathematical difficulty remains unchanged 
and has to do with the strong non-locality of the dynamics. It is built in in 
the definition of the sandpile process and it seems to be essential for some 
of its most interesting aspects. The effect of adding a sand particle at one 
site can affect the whole system. This implies that the standard machinery 
to construct a Feller process (Hille-Yoshida theorem, martingale problem) 
does not work. Similar problems were encountered in the study of long range 
exclusion processes, see |J, ||, @. Here the situation is even more exciting 
as will be seen from the convergence to the stationary measure in a finite 
time and from the relation between the process and its formal generator. 
This relation is only 'pointwise' for a certain class of 'decent' configurations 
(in particular the process does not have a generator in the standard sense). 
One finds analogies with some recent work on weakly and almost Gibbsian 
measures where the interaction energy is only absolutely summable on a full 
measure set of configurations but not uniformly, see jLl], [12| . 

Our paper is organized as follows. In Section 2 we introduce the model 
and the formal generator of the sandpile process. In Section 3 we prove how 
to construct a Markov semigroup, and thus a Markov process which corre- 
sponds to this formal generator. The main technical tool in this construction 
is monotonicity. In Section 4 we show the precise relation between the semi- 
group and the formal generator. In Section 5 we study the invariant measure 
of the process and show that every initial configuration converges in a finite 
time to this invariant measure. 



2 Definitions. 

2.1 Configurations. 

The state space of the one- dimensional sandpile process is Q = {1,2}^. 
Elements of Q are called configurations and are denoted by rj, £, (. The state 
space is a compact metric space in the product topology. The value r)(x) = 2 



3 



is called the threshold value. Sites for which r)(x) = 2 are called critical sites 
in the configuration rj. The set of all finite subsets of TL is denoted by S. We 
define Qf to be the set of configurations with a finite number of critical sites, 
i.e. 

n f := {q G n : ^ 1 ({2}) G S}. (2.1) 

There is a one-to-one correspondence between S and Qf via A — > t]a, where 
rjA denotes the configuration whose set of critical sites is A, i.e. ^ 1 ({2}) = A. 
The configurations in flf will play an important role in the construction of 
the process (cf. below). A function / : Q — > 1R is called local if it depends 
only on a finite number of coordinates. Every local function is continuous 
and every continuous function is a uniform limit of local functions. 

2.2 Toppling transformation T{. 

Given a configuration rj G and a lattice site i G 2, Tj?? represents the 
configuration obtained by adding one sand particle at site i and letting the 
system "topple" until a (stable) configuration (i.e. an element of Vtf) is 
obtained. This toppling is described in words as follows: If rj(i) = 2 and 
one sand particle is added at site i, then from that site two particles will be 
removed, one to the left neighbour % — \ and one to the right neighbour i + 
If say r](i + 1) was already 2, then site i + 1 will also "topple", giving one 
particle to i, and one to i + 2. This goes on until a configuration in which 
v(j) < 2 for all j G TL is obtained. This succession of topplings at different 
sites caused by adding one particle to a critical site is called an avalanche, 
and thus Tit] represents the configuration after the avalanche. If r](i) = 1, 
then Tii] differs from r] at site i only (TiT)(i) = 2), but if r)(i) = 2, then 
differs from rj at three lattice sites which depend on i and rj in a non-local way. 
This corresponds exactly to the dynamics of the standard abelian sandpile 
process in a finite lattice region. 
More precisely we introduce 

A;+(*,?7):=inf{j>0:7^+j) = l}, (2-2) 

and 

k-{i,ri) := inf{j > : V (i - j) = 1}, (2.3) 

where inf := +oo. We denote by e, the mapping from TL to {0, 1} which is 
one at site % and zero at all other site. We distinguish five cases: 

1. k + (i,i]) = 0, i.e. r)(i) = 1, then 

= + e t . (2.4) 
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2. k + {i, 77) > 0, k + (i, 77) V k (i, 77) < 00, then 

Tifj — T] + ej_|_/ c +(j j ^) + &i—k-{i,r}) ~ e i+k+(i,ri)-k~(i,ri) (2-5) 

In what follows we abbreviate 

F ■= i + k + (i,r]) - k~(i,rf). (2.6) 

3. k + (i,r]) = oo, k~(i,rj) < oo, then 

T,r] = r] + e^k-^) (2.7) 

4. k + (i,rj) < oo, k~(i,rj) = oo, then 

TiT] = r] + e i+k +^ (2.8) 

5. k + (i,rj) = k~(i,r)) = oo (i.e. 77 = 2), then 

= 77. (2.9) 

Cases 3,4,5 above are just conventions to define Tiq for all 77 e f2, and to 
write down the formal generator, but they are not relevant in the construction 
of the sandpile process. That these choices are 'correct' will be seen later 
once it is shown that the configuration 77 = 2 is invariant for the constructed 
process (cf. Proposition |5.1|). 

The explicit and simple expression of Tit] is typical for the model in di- 
mension one. In higher dimensions TjT] will become a much more complicated 
function of rj. Note that the transformations Tj and Tj commute, i.e. 

7^(77)) = TXTifa)), Vz, j ei,y v e n. (2.10) 

For this reason the sandpile model has been called abelian in the physics 
literature, cf. 0. 

2.3 Monotonicity. 

In the construction of the sandpile process, monotonicity will play a crucial 
role, just as in the construction of the long-range exclusion process, cf. 0, 
||. For more details on monotonicity, we refer to ||10|| , Chapter 2, Section 2. 

For 77, £ G Q we define 77 < £ if i](x) < £(x) for all x € 7L. A function 
/ : Q — > IR is called monotone if 77 < £ implies f(rj) < /(£), for all 77, £ G 
fi. We denote by the class of all monotone bounded Borel measurable 
functions. A Markov process {f] t : t > 0} on Q with Markov semigroup 
{S(t) : t > 0} is called monotone if / G M. implies S(t)f G M.. Usually one 
proves monotonicity of a Markov process by explicitly constructing a coupling 
of path space measures such that for rj < £ one has P v ^(r] t < £ t , Wt > 0) = 1. 
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2.4 Formal generator. 

In sandpile processes, the avalanches are immediate. Once the vertex i is 
chosen, the transformation r\ — > happens instantaneously. This is what 
one calls infinite separation of time scales in the physics literature, and it 
has been argued as essential to have a so-called self-organized critical state. 
The formal generator of the sandpile process has therefore to be written 

as 

mv) = Y,mv)-f(v)]- (2.H) 

iel 

Since L is of course not the generator of a Feller semigroup, we shall have 
to give a precise meaning to this description of the dynamics, and to clarify 
how to "associate" a process to L. 

We can rewrite the formal generator 

Lf(r,)=L b f(r,)+L a f(r l ), (2.12) 

where 

Lbfiv) = E X(v(i) = W(V + <*i) - f(v)} is the "birth part", (2.13) 
iel 

L a f(v) = E Xfa(i) = 2 )[/(^77) - f(r,)] is the "avalanche part", (2.14) 
iel 

and x denotes the indicator function. The operator L b corresponds to a pure 
birth process (independent births with rate one at each lattice site) and thus 
is a well-defined generator of a Feller semigroup. The avalanche part however 
gives technical problems related to the nonlocality of the transformation T{. 
One cannot run the standard Hille-Yoshida program, which would give a 
Feller semigroup associated to L b + L a . Indeed, the process we shall associate 
to L will be constructed by monotonicity and turn out not to be Feller. 

Remark. The processes we consider here run in continuous time while 
it is a standard procedure to describe the finite volume sandpile process in 
discrete time. However, since both updating mechanisms are sequential, the 
difference is unimportant. 

3 Construction of the sandpile process. 

We proceed in three steps: 

1. Definition of a process associated to L a (only avalanches) on Qf. Con- 
struction of a coupling showing that this process is monotone. 



6 



2. Definition of a process associated to L a + births in the finite interval 
[— n, n] on Qf. Construction of a coupling showing that this process 
is still monotone. This will give us a semigroup defined on / G M. 
by S n (t)f(r)) = lim^/ e n yj S n (t)f(r)'), where liny^ denotes the limit 
along sequences {rj' n ,n G IN} C Qf such that i]' n < rf n+1 . 

3. Monotonicity of the semigroups S n (t) in n, i.e., for all / G Ai, rj G Q, 

t > 0, S n+1 (t)f(r}) > S n (t)f(rj). 

We shall finally define the sandpile process as the process associated to the 
semigroup 

S(t)f(rj) = lim lim S n (t)f(rj'). (3.15) 

3.1 Step 1: Avalanche part of the generator. 

For i] G flf, A(rj) := rj~ l ({2}) G S and we can construct the process with 
generator L a from the non-exploding continuous time Markov chain {A t : 
t > 0} on S with generator 

LJ(A) := J2 {f(A \ {<"} U {t + fc + (7, 77)} U {i - k~(i, 77)}) - f(A) , 

ieA 

(3-16) 

for / : S — > H bounded. Indeed, IE A |A t | = |A|e* (where |A| denotes the 
cardinality of A) so that A t G 5, for all t > 0, P A -a.s. The process with 
generator L a (avalanche part) is thus defined on Qf by ^ := 77^, where 
{At '■ t > 0} starts from A = i]^ 1 ({2}). We show that this process is 
monotone via the coupling of rjt, ( for 77, £ G f2/ with £ < 77) described as 
follows. When a sand particle is added in 77 at site i G 2 where 77(2) = 2, 
an avalanche creates a 1 at site i v = i + k + (i,rj) — k~(i,rj). Define now 
0(7, 77, as the unique site on Z such that a sand particle added at that site 
in configuration £ would create a 1 in £ at the same site i v . More precisely, 
if £(7") = 2, 

0(7,77,0 = ^ + ^(^,0-^(^0, (3-17) 
and if £(i' J ) = 1, we define 0(7,77, = 00 and = £. We now write the 
coupling generator on Qf x flf for initial configurations (77,0 with 77 > £ as 

L c a f(v,0 = E X(v(i) = 2)[/(T i 77,T^ ) - /fa, 01- (3-18) 

Proposition 3.1 1. On Vtf x f2j the generator L c a defines a Markov pro- 
cess. For 77 > i/ie pai/i space measure P v ^ is a coupling of P v and 
P^ (the path space measures of the processes with generator L a starting 
respectively from 77 and O- 
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2. Forr]>i we have P v ^(r] t > Vt > 0) = 1. 



3. The process with generator L a is monotone on Vtf. For the semigroup 
S a {t) := exp(tL a ) associated to L a , f G M. implies S a (t)f G .M. 

Proof: To check that the process with generator L c a has the right second 
marginal for rj > £, remark that for every j G 7L such that = 2, there 
exists a unique i G 7L such that 0(i,r/,£) = j (L„ does not define a coupling 
of P v and for general couples (77, £) but only for the ordered couples, i.e. 
(77, £) with 77 > £ and 77, £ G fi/). For points 2 and 3, note that by construc- 
tion any transition in the coupled process preserves the order. ■ 



3.2 Step 2: Avalanche + birth in a finite interval. 

We now include the birth part of the formal generator L, only on the interval 
[—n, n}. We have to prove that the associated process on flf is still monotone, 
since births and avalanches do not commute. For this, the generator 

L (n) f(v) = Ex(l(i)=2)[/(rffl)-/(?)] 

n 

+ £ x(r](i) = W(Ti7l)-f(ri)} (3.19) 

i=—n 

defines a continuous time Markov chain on S, and thus a process on Qf. We 
now couple two versions of this process starting from 77, £ G f2/ such that 
V > £, with the generator 

L c > {n) f(v,H) = ExW0 = W)=2)[/(3li?,W)-/('J ) O] 
+ £ = 2)x(£(0 = 1) [f(v, m - f(ri, 0] 

n,n] 

+ e xWi) = i)x(ew = i) [mv, m - nv, o] 

n,n] 

+ Ex(»7(0 = 2)x(e(«" ,| ) = i)[/(r i i7, 01- (3-20) 

i 

This coupling is described as follows: To each site i G TL we associate two 
independent Poisson clocks, C and C. C indicates the "avalanche event- 
times" and C the "birth event-times". When C rings at % and r\{%) = 2, 
we simultaneously add a sand particle at % in 77 and at 77, £) in £ (if 
(frihViO = °°5 we have no transition in £). When C rings at 2 and r/(i) = 1, 
then, since 77 > £, = 1 and we add a sand particle in both £ and 77 
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simultaneously. If however rj(i) = 2 and = 1, then we only add in £ at 
site %. Like in Section 3.1, this construction implies 

Proposition 3.2 1. The generator L c '^ defines a Markov process on 
Qf x Qf. For r] > £ its path space measure P^ v '^ is a coupling of 
P% and Pj, the path space measures of the processes with generator 
starting respectively from 77 and £. 

2. Ifr]> f, then P^ ] {r] t > Vt > 0) = 1. 

3. The process with generator is monotone onflf. For the associated 
semigroup S n (t) := exp(tL^), / e M. implies S n (t)f £ M.. 

3.3 Step 3: Avalanche + births in different intervals. 

This final step is the easiest: By allowing births on a larger interval, we can 
only obtain larger configurations. We couple the processes with generators 
]Jn+i) anc j jjn) through the following generator, where 77, £ e Qf, 77 > £. 

+ J: xm=2)xm = w(v,m-f(v,o\ 
+ 2 x(i7(o = i)x(e(o = i)[/(^,TiO-/(i7,o] 

ie[— n,ra] 

+ £ X fa(i) = i) x (e(<) = i)[/(^,0 - 7(77,01 

71— i,n+i} 

+ Y,xW) = 2 MW) = 1 MTv,Z)-f(v,Z)]- (3.21) 

Proposition 3.3 1. The generator L c ' n+1 ' n defines a Markov process on 
Qf xQj. For i] > £ the path space measure Pn+i n ^ s a coupling of P^ +1 
and Pi*. 

2. For V >t P^ lin ( Vt > Vf > 0) = 1. 

3. For f EM, riE n f , we have S n+1 (t)f(r]) > S n (t)f(r]). 

Proof: By construction any transition in the process with generator L c,n+1,n 
preserves the order. To see point 3, for / G M. 

S n+1 (t)f( V ) = Ej^/fo) = E^ 1 , n /(7 7t 1 ) > niij^t) 

= WPJ{ m ) = S n {t)f{r 1 ). (3.22) 
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3.4 Construction of the process. 

We can now proceed to the definition of the sandpile process as a well-defined 
Markov process with formal generator L. For functions / G Ai, on Q we 
define the semigroup associated to L n by 

S n (t)f(v) = , lim S n (t)f( V >). (3.23) 

In ( |3.23| ) the limit is taken along the increasing sequence i]' n defined by 
T]' n (x) = r)(x) for \x\ < n, and ^(x) = 1 otherwise, but any other increasing 
sequence rf n f 77 gives the same limit. Then, for / G .M and r/ e O, 

S n+ i{t)m > S n (t)f( V ). (3.24) 

Since S n (t)f(r}) < ||/||oo (where [|/||oo = sup^ en |/(^)|), for / G and 
G we can define 

S(t)f(r,) = lim S n (t)f(ri). (3.25) 

n|oo 



Because each ^(t) is a semigroup on functions on Qf, we have, using ( p.24|) , 
that S(t)l = 1 and that S(t)f > if / > 0. Moreover, for f E M, 



S(t + s)f = \imS n (t + s)f = \imS n (t)S n (s)f 

njoo n\oo 

< lim S n (t)(S(s)f) = S(t)S(s)f, (3.26) 

n| 00 

and for all m G IN 

S(t + s)/ = lim^^^^/^liminf^^)^^)/) 

njoo njoo 

> S m (t){limMS n {s)f) = S m {t)S(s)f, (3.27) 

njoo 

so by taking the limit as m j 00 we derive the semigroup property for S(t). 
This semigroup defines the finite dimensional distributions of a Markov pro- 
cess; by Kolmogorov's theorem, there is a unique Markov process {r]t : t > 0} 
such that S(t)f(j]) = JEPf(j]t). We call this process the one- dimensional 
sandpile process (SP). We denote its path space measure (i.e. a measure on 
starting from 77 by -Pyp. 

Remarks. 

1. We cannot conclude that Pg P has a version concentrating on D([0, 00), Q) 
(the set of cadlag paths on Q) for all initial 77 (cf. Corollary 5.1 below). 
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This question is related to the following one: For which r] G do we 
have 

Bm£(f)/fa)=/fo) (3.28) 

for all local /? It is not even clear that this limit exists. In contrast, 
this is automatic for long range exclusion: But in our case, because the 
values of a configuration in a finite volume are influenced by transitions 
outside this volume, the estimate (2.3) of U does not hold. 

2. We cannot expect that the formal generator L is "really" the generator 
of the process. However we do prove in Theorem 4.1 that for some class 
of "decent" configurations and local / 

]im±(S(t)f( V )-f(ri)) = Lf( v ). (3.29) 

This situation of pointwise convergence for a set of "decent" config- 
urations, rather than uniform convergence, can be compared to the 
situation of "weakly Gibbsian measures" where the potential is only 
convergent on a set of good configurations, cf. |TTJ,Jl2|. In the fol- 



lowing section we will deal with ( |3.29j ) and show that for short times 
(depending on 77) and for "decent" configurations (a set of measure 
one for every ergodic probability measure on Q, except 62 , the measure 
concentrating on 77 = 2) we actually have 

S(t)f(v) = E -<L n f{v), (3-30) 



which is of course much stronger than ( |3.29| ) (but note that we do have 
Q3.30D for the pure birth process). 



In order to find an invariant measure for the sandpile process, we cannot 
proceed in a standard way. Since the process is not Feller (see Corollary 
5.1 below), a weak limit point of ^ Jq /j,S(t)dt does not need to be 
invariant, and / Lfdfi = does not imply that fi is invariant. In 
Section 5 we shall show that the only invariant measure is 82 and that 
every initial measure converges to 82 in a finite time. This explosive 
convergence to the invariant measure is related to the non-locality of 
the dynamics and to the fact that we continue to add sand at rate one 
at each lattice site (cf. @). 
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4 Decent Configurations. 



As we mentioned in Section 2, the transformation Tj is non-local. Neverthe- 
less we still have some kind of locality: Let / : il — > 1R be a local function 
depending on coordinates in A E S. Put A~ := min(A), A + := max(A), 
then we have 

TJiv) ~ f(v) = 0,\ftEl\[A-- k-(A-, V ), A + + k + (A + , 77)]. (4.31) 

This means that the effect of adding a sand particle is only felt in the region 
bounded by the first one to the left and the first one to the right of the 
dependence set of /. Loosely speaking, Tj/ has a dependence set which is 
finite but dependent on the configuration. 

We now formalize these statements. Let Qi denote the set of configura- 
tions with an infinite number of ones to the right and to the left of the origin, 
i.e., 

Q l := { V G n : E(2 - v(i)) = E( 2 " = +~}- ( 4 -32) 

;<o i>o 

For t) G Hi we order the sites % G 7L for which r)(i) = 1 as follows 

^^({l}) := {X^rj) :ieTL} (4.33) 

where X Q (r]) := min{i > : r}(i) = 1} and i < j implies Xiijf) < Xj(rj). 
Then we define 

/o =(x_ 1 ,x ]nz 
h =(x ,x 1 ]nz 

l_ x = (X_ 2 , X_i] n 1, etc. (4.34) 

In this way we view configuration rj as a sequence of intervals, where each 
interval is bounded by two sites i,j G 7L, i < j, for which rj(i) = r)(j) = 1. A 
function / : Q 1 — > IR will be called N— local around k G Z if 

/fa) = /^(0, ^ U^^^j • ( 4 - 35 ) 

We say that / is A-local if it is A-local around the origin. 
Examples. 

1. f(rj) = rf(i) is 0— local around i, and every local function is A-local for 
some N . 
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2. f(j]) = |/fc(^)| is A;— local for every k, but is not local. 

3. f{rj) = Sj; 6 2 e ~ v( x )j while continuous, is not iV-local for any N. 
Proposition 4.1 Let f : Qi — ► 1R be N— local then 

N+l 

f{T iV )- f{ V ) = for alii e 1\ [j 1,(77), (4-36) 

j=-N-l 

and Lf is (N+l ) -local. 

Proof: If i G ujL + _Vi J j(^)> then ; for / AMocal, f(T^) is (N + l)-local. 
Otherwise Tit] — r\. ■ 

Proposition 4.2 Let f : Qi — > 1R fre N— local and bounded, then we have 

\L n f( V )\ < (\I-n-u(v)\ + ■■■ + |W(^)|) n 2 n ||/||oc- (4.37) 

Proof: We put N = for the sake of simplicity, the case N > is treated 
in the same way. We proceed by induction. For n = 1 

Lf( V ) = Yllf( T iV)-f(v)], (4-38) 



and, by Proposition |4~l| , this sum runs over i G I-i(t}) U Ioiv) U A(^), A/ is 
1— local, and 

\Lf( V )\ < (IM^I + 1/0(17)! + 1/1(17)1)211/1100. (4.39) 
Suppose now that for n > 1, L n / is n— local and 

|/ n /(^)l < (|I-«fa)| + • • • + \In(v)\n n \\f\\oc. (4.40) 

Then we have 

L n+1 f(v) = Y\L n f{T iV ) ~ L n f( V )], (4.41) 

i 

where the sum runs over i G I- n -\{rj) U . . . U I n+ i(r)), since L n f is n— local. 
Therefore, for such an i, 

in n || f\\ 

00 



|/"/(^)l < [\I-n(li V )\ + . . . + \I n (T iV )\] n T 

< [1/^(77)1 + ... + |/„ +1 (77)ir2l/|| 0O . (4.42) 

Thus 

\L n+1 f(v)\ < [\I- n -x(v)\ + --- + (v)\) n+1 2 n+1 \\f\U (4.43) 

■ 

The following lemma is a direct consequence of Cauchy's formula for the 
radius of convergence of a power series. 
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Lemma 4.1 Let {a n : n > 0} be a sequence of positive real numbers such 
that limsup^^ a n /n = a < oo. Then the series J2^=ot na n/ n ^ converges for 
\t\ < 

I I a.e 

We now define the set of "decent" configurations: 

Q dec ■- { v e Q l : a( V ) := limsup l J -"^)l + • • • + l J "(^)l < ^ (4 44) 

Observe that for rj G f2/, if |z| is large enough, = 1, hence a{ri) = 1. 

Therefore fi/ C ^ci ec - More generally, if the density of ones 

pf v ) = l im ^=-n(2-r/(i)) ( } 

HK " «-oo 2n+l v ' 

exists, then a(r]) = l/p(rj). This implies that for every probability measure 
\i ergodic under spatial translations and for which /i(2 — 77(C))) > 0, the set 
f2^ ec has //-measure 1. It is also clear that ^dec * s a translation invariant 
set in the tail field, i.e., if i] G ^oi ec and £ differs from rj in a finite number 
of lattice sites, then ( G ^dec 

Theorem 4.1 Let f be N— local for some N G IN and rj G ^d ec , then for 
t < 4et ^ £/ie series Y,T=o * L n f ^ converges absolutely and equals S(t)f(r)). 
In particular, L is the "pointwise generator" of the semigroup, i.e. 

^ mm-m = Lm ( , 46) 

Proof: The absolute convergence of the series follows from Proposition 4.2, 
Lemma 4.1 and ( 4.44[) . By definition, for an iV-local / G A4 

S(t)f( V ) = lim lim S n (t)f( V '). (4.47) 

n\ oo rf &lf ,r) \rj 

For rf G Qj, a(r)') = 1 and 

oo 

S„(f)/(i/) = £*V* ! (A0W)- (4.48) 

/c=0 

Moreover, for every /c G IN, 77 G f2: 

hm lim (L n ) fe /(V) = L k f{rj). (4.49) 

njoo 77 '&lf,r) Jr) 

We can bring in the limits lim^oo lim^/^ into the sum of ( 4.48|) by the dom- 
inated convergence theorem together with Lemma 4.1 and the inequalities 

< (\I- N - k (v')\ + ... + \I N+k (v')\) k 2 k \\f\\oo 

< (\I^ k (r ] ) + ... + \I N+k (7 1 )\) k 2 k \\f\U (4.50) 
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for 7]' < T], T]' G Of. 



Remark. Theorem |4.1| shows that the infinite volume sandpile process 
is indeed a natural extension of the finite volume standard abelian sandpile 
model. More precisely, if P n denotes the transition probability matrix of the 
discrete-time sandpile model in the finite volume A n = [— n, n] D Z, then we 
have for 77 G fi^ec and t small enough (depending on 17), 

limPW/^^W/^, (4.51) 

n\oo 

where rf 1 denotes the restriction of rj to A n . 



5 Invariant measure for the sandpile process. 

In this section we show that after a finite time, the process reaches 62 , which is 
therefore its only invariant measure. Nevertheless we first prove Proposition 
5.1 to illustrate some "pathological" behavior of the process (cf. Corollary 
5.1). 

Proposition 5.1 62 is invariant for the sandpile process. 

Proof: Denote by 2 the configuration rj = 2. Notice that this is not a decent 
configuration. Thus we have to prove that 

S(t)f(2) = /(2). (5.52) 

Without loss of generality we put f(r)) = n{0) - 1. Clearly, S{t)f(2) < f(2), 
so we are left to show the opposite inequality. By definition ( |3.25| ), 



S(t)f(2) = lim lim S n (t)fW) 

n^oo n' d\lf,ri'\ri 

> lim S a (t)f( V ') 

•q' &lf ,rj'\rj 

= lim P An (0 G A(t)). (5.53) 

nfoo 

Here P An denotes the Markov measure on (S' ' 00 - 1 associated to the process 
with generator L a (only avalanches), starting from the set A n := [— n, n] H Z. 
The set A(t) will always be of the form I(t) \ {x} where I(t) is an interval 
which grows in time. The probability that {x} = {0} will thus be of the 
order (\A n \ + /(t))" 1 where f(t) is some increasing: function of t. Let us 
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make this precise. Denote by Ti, r 2 , . . . , r n the successive event times of the 
Markov chain {A(t) : £ > 0}. Then we have 



1 



14 r 

P A "(0^(r 2 )) = E = 5 ) = nTTXT' ( 5 - 54 ) 



and analogously 



P A "(0^A(r n ))^ |At|+ 1 w _ 1 - (5-55) 



Hence we conclude that for alH > 



Therefore 



limsupP A "(0 ^ A(t)) = (5.57) 

n^oo 



which proves the claim. 



Corollary 5.1 The sandpile process is not Feller, and for some initial con- 
figurations the process has no right- continuous version. 

Proof: As in the preceding proof, we have that for all rj e £1 such that 
?7 _1 ({1}) G S (finite number of ones) 

P V sp(Vt(x) ^ 2) = 0, Vt > 0. (5.58) 

Therefore there is no right-continuous version of the process starting from 
these configurations (cf. Remark 1 of Section 3.4). ■ 
The following theorem will be proven at the end of this section. 

Theorem 5.1 In the sandpile process there exists a finite T > such that 
for all t > T , i] E Q, 

PspIvM = 1] = 0. (5.59) 
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This result seems intuitive: Since we continue adding sand at each lattice 
site with rate 1, after a time of order unity, every site will have received an ex- 
tra sand particle. Once the configuration is full, it remains full. Nevertheless, 
this intuitive picture is somewhat dangerous, since we are working in infinite 
volume, and it could be possible that sand particles "disappear" at infinity. 
Indeed, the same statement is not true anymore in higher dimensions: Al- 
though we still expect convergence to the invariant measure in a finite time, 
this invariant measure will not be the Dirac measure concentrating on the 
"full configuration" . 

5.1 The finite volume sandpile process. 

For the finite volume sandpile process in [—n, n] , we only add sand particles 
in [— n, n], and at the boundaries we put fixed ones (i.e., sand falls off the 
boundary and disappears when the boundary site topples). More formally, 
for rj G Q, we denote by rf 1 the configuration which is equal to r\ in [—n, n] 
and is identically 1 outside that interval. Define 

T n , V := [T tV n } n , (5.60) 

fy,n := {n e n : 77O) = 1, Vj G Z \ [-ra, n]}. (5.61) 

Note that Qf — U n Q/ jn . The finite volume sandpile process (FVSP) is the 
pure jump process on n with generator 

n 

£nf(v) = E U( T n,v) - f(v)\ (5-62) 

i=— n 

5.2 Poisson representation. 

Let {N t : t > 0} := {NJ:,i G 7L, t > 0} be a collection of independent rate 
one Poisson processes. For S G 1N we define 9 n {E) to be the configuration 
obtained from "toppling" in [— n, n] and fixing ones outside, i.e., 

1 (5.63) 

where i denotes the configuration rj = 1, and T^" 1 means E(i) — 1 iterations 
of T n>i . 

Proposition 5.2 The process {9 n (N t + rj) : t > 0} is a version of the FVSP 
starting from r\ G Vt^ n (N t + rj means coordinatewise addition). 



*n(3) = 



2(0-1 
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Proof: Notice that the transformations T n i commute for different i. There- 
fore, first adding all the particles and then letting the system topple is equiv- 
alent to adding the particles one by one and letting the system topple each 
time. More formally, using ( |5.63| ), 

(±Ef(9 n (N t +v))^j = C n f( V ) (5.64) 
for r] G % n . ■ 



5.3 Invariant measure for the FVSP. 

Since the FVSP is a pure jump process on the finite set fi/, n , it converges 
exponentially fast to its unique stationary measure /j n . 

= T^l E S « ( 5 - 65 ) 
is the uniform measure on the set TZ n of recurrent configurations: 

n 

n n = {r 1 e Q f , n : EC 2 " Vii)) < 1}, (5-66) 

i.e., those configurations in Qf >n containing at most one site i G [— n, n] such 
that = 1. For all i G [— n, n] fl Z, ^ n T n ,i = Hn- More precisely for every 
pair (//, C) £ Tin x 7l n , V C there exists a unique lattice site i G [— n, n] 
such that T n i ri = (. Therefore, 

r 1 n 

(C n f)gdfi n = EiSlE [f(T n , lV )-f( V )]g( V ) 

V£Kn I ™l i=-n 

= T^i E 1/(0 - /fo)Ufo) 



f(C n g)dfj, n . (5.67) 
so that /!„ is actually reversible for C n . 

5.4 Convergence to equilibrium. 

The key inequality in the proof of Theorem 5.1 is: There exists T > such 
that for all t > T 

P{6 n {N t + 77) (0) = 1) < Ce- nX + C'/n. (5.68) 
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Indeed, after a finite time T in the Poisson process we have for large n 

n 

n t~ 2nT - ( 5 - 69 ) 

i=— n 

For T large enough this represents a serious "excess" at most lattice sites in 
[— n, n\. If we let this configuration topple, then we will obtain a recurrent 
configuration, as we show in the following lemma. 

HI 

Lemma 5.1 Let rj G IN be such that 

Ln/2j 

E ^)>12n, (5.70) 

i=-\n/2\ 

then n (r)) G lZ n . 

Proof: First of all note that we put the "excess" of particles sufficiently far 
from — n and n, so that not too much sand will be lost at the boundaries. We 
call a site i G [—n, n] n 7L an excess site if r)(i) > 2. For the sake of simplicity 
we suppose n even. Put 

m = min{i G [-n/2, n/2] n Z : r)(i) > 2} 

M = max{i G [-n/2, n/2] n Z : 77(f) > 2}. (5.71) 

Because of the monotonicity property 

77, C e IN Z , C > V, and 9 n (n) G ^ n imply 6 n {() G (5.72) 

we can replace our configuration i] by rj' defined by: f]'(i) = r](i) for % G 
[—n/2, n/2] fl TL and n'(i) = 1 otherwise. We now let the system topple in 
the following way: All sites in (m, M) topple and the sand arriving at m or 
M will be stocked at these sites (i.e. sites m and M do not topple). We then 
end up with a configuration rj" with two excess sites at m and M and since 
no sand is lost in this toppling: 

n"{m) + n"(M) > lOn. (5.73) 

We change again the configuration only retaining the site with the maximum 
excess, say M, and replacing all the rest by l's, this configuration is denoted 
by rj'". The excess of r/'" at M is not less than 5n. So we are left to prove 
that a configuration with one excess site of this type topples to a recurrent 
configuration. The case M = is trivial. Suppose M > 0. After n — M 
topplings we get a configuration rf v with an excess site at M (an excess not 
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less than 3n), rf v (i) = 2 for all i G [2M — n,n], n lv (i) = lmbox for all 
i 6 2 \ [2M — n, n] . From that moment on, sand particles are lost at the 
boundary by further toppling. It is easy to see that the worst situation is 
when M = n/2. In that case, define the configuration rf by rj v (i) = 2 for 
i G [0, n] and rf{i) = 1 otherwise; since the excess at site i = n/2 in rf v is 
not less than 3n, it suffices to show that {T n>n i^f n {rf) £ Ttn- 

By adding sand particles at site n/2 in rj v , we shall always obtain a 
configuration in the set T n = U_ n < a <or n ,a, with 

n 

r„, a = {ne Q f , n : n(i) = 1 V? < a, and ]T(2 - < 1}. (5.74) 

i=a 

That is, for n G r njQ , there exists a "leftmost" site a < for which 77(0) =2 
and the interval [a, n] contains at most one site i for which rj(i) = 1. When a 
is the position of the leftmost two in the configuration rj, then by adding three 
sand particles at site n/2 and letting the system topple, the position of the 
leftmost two becomes a — 1. I.e., for a < 0, 77 G r„ ia implies (T„ in / 2 ) 3 G r n>a _i. 
By this last property, since T n _ n = TZ n , {T ntn/2 ) 3n n v G TZ n . ■ 

The following lemma is a standard large deviation estimate for sums of 
independent random variables. 

Lemma 5.2 For every strictly positive t ande, there exists a constant A(e) > 
such that 

( \n/2\ \ 

P\ N?<n(t-e) <exp(-nA(e)). (5.75) 

\x — l»/2J / 

Proof of Theorem 5.1: We show that the statement is true for T = 14. 
By a coupling argument as in Section 3, we get that the FVSP 9 n (N t + rj) is 
stochastically dominated by the sandpile process rjt (since we put fixed ones 
at the boundaries of the interval [— n, n] in the FVSP). Put < 5 < 1. By 



Lemmas |5.2| and |5.1| we know that for t > 14, there exists A > such that 

P(e n (N t - S ) & K n ) < Cexp (-An). (5.76) 

Hence we obtain the inequalities 

PMvt(0) = l) < P(6 n (N t + n)(0) = l) 

< P[9 n (N t + V )(0) = l\9 n (N t _ 6 + r ] )eTZ n ] 

+ Cexp(-Xn). (5.77) 
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Conditioned on {( := 9(N t _s + v) e ^-n}, the distribution of 9 n (N t + rf) is 
uniform on a set of the form lZ n \ {£'} if there has been at least one Poisson 
event in (£ — 5, t). Therefore, since there is only one rj G TZ n with 7/(0) = 1, 



P [6 n {N t + ri){0) = l\6 n {N t - S + ri)ieK n ] 



-<5\2n+l 



2n + 1 

Combining (|5.78 ) with ( 5.77|) we arrive at 

1 



P2pfo(P) = i] < 



-5\2n+l 



2n + 1 v 

which concludes the proof of the theorem. 



Ce 



-nX 



(5.78) 



(5.79) 
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